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Abstract 

The impact of the recently evaluated 'irreducible' contributions to the resummed next-to-leading 
logarithmic small-a; anomalous dimension 7 99 is evaluated for the unpolarized parton densities 
and structure functions of the nucleon. These new terms diminish the gluon distribution and are 
found to overcompensate the enhancement caused by the resummed leading logarithmic small-x 
anomalous dimension and the quarkonic contributions beyond next-to-leading order. 
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The structure functions of the proton in unpolarized deep-inelastic scattering (DIS) show a steep 
rise towards very small Bjorken-x values, which becomes stronger with increasing resolution Q 2 . 
This behavior is a consequence both of the shape of the non-perturbative quark singlet and 
gluon initial distributions - xS and xg rise roughly like x~ ' 2 for a starting scale Qq ~ 4 GeV 2 
- and of the form of the evolution kernels governing the renormalization group equations of the 
mass singularities. The anomalous dimensions 7^ for the evolution of the parton densities, as 
well as the Wilson coefficients C n of the structure functions, contain large logarithmic small- 2 
contributions. In order to arrive at a reliable theoretical framework at very small x, the resum- 
mation of these terms may be necessary to all orders in the strong coupling a s (Q 2 ). For the 
unpolarized singlet case considered here the dominant contributions take, in Mellin-iV space, the 
form a l s +k /(N - l) k , I > 0. 

The resummation of 7 S9 and 7 OT in the leading small- x approximation (Lx), I = 0, was 
performed long ago 7^(iV) is obtained as the solution of 

l=j£lXb(lS>) (1) 
with a s = C A a s (Q 2 )/n, C A = N c = 3, C F = 4/3 and 

Xo(T) = 2V(1) -V<7)- V>(l-7), (2) 

furthermore 7^ = (Cf/Ca) Jg°J ■ The quark anomalous dimensions j qg and 7 g9 , on the other 
hand, receive contributions for / > 1 only. The I = 1 terms were derived in ref. 0, together 
with the corresponding resummations for the coefficient functions Ci and C^. The large effects 
of these quantities on the small- x behavior of the DIS structure functions were subsequently 
studied in detail [3-7]. In those investigations the small- x resummation of the gluon anomalous 
dimension ^ gg to next-to-leading order small-x (NLx) accuracy could not be taken into account. 
This resummation has now been performed for the quarkonic contributions [8-11] proportional 
to the number of quark flavors Nf. Recently also the 'irreducible' gluonic terms oc Ca have been 
derived |fl2| , i.e., those contributions which are energy-scale independent in the framework of 



ref. |10[ underlying that calculation. The corresponding terms of 7 9g , however, still remain to be 
determined. 

In this note we investigate the impact of these new resummed contributions to 7 99 on the 
evolution of the parton densities and the proton structure functions F 2 (x,Q 2 ) and F L (x,Q 2 ), 
for the first time including calculated subleading terms into the renormalization group analysis. 
Hence the comparison of the results to the findings of previous studies [3-7] should allow for 
improved estimates of the convergence of the small-x resummation approximation, despite a 
fully quantitative, scheme-independent NLx analysis not being possible at present. 

As will be demonstrated below, the effect of the new contributions to the resummed anoma- 
lous dimension 7 99 is very large and opposite to that of the previously known resummed terms. 
This implies, already at the present stage, considerable changes particularly for gluon-dominated 
quantities, which partly modify conclusions obtained in previous numerical investigations [3-6]. 
A detailed account of the solution of the evolution equations in the presence of all-order anoma- 
lous dimensions and coefficient functions will be given in a forthcoming publication [|T^ . 

As shown in ref. [p| the larger eigenvalue of the singlet anomalous dimension matrix, 7+ (TV), 



may be obtained in the Qo scheme |14| as the solution of 



[Xo(7+) + a.Xi(7+)] , ( 3 ) 
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where the second term is the sum of 
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The function ^(7) in eq. (H) is given by 
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with ai = 0.72, a2 = 0.28 and 03 = 0.16 |12| . From these results the irreducible NLx-contribution 
to 7 99 is then inferred by 
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Our subsequent numerical analysis will be performed in the DIS factorization scheme. Here 



7^ is represented as 



lgg,T>lS 
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with R(a s ) defined in ref. p[ and /5 = (11/3) C^ — (2/3) iVy Tables of the expansion coefficients 
dgg'j^' b \ dgg k and ffc may be found in ref. |L3[]. Here we list for brevity only the numerical values 



of the first 15 coefficients b k and bj! A1 ' for the Lx and NLx series for Nf = 4, see Table 1. Note 
that the new terms bfj'^ agree with the corresponding results from fixed-order perturbation 
theory already taking into account the 'irreducible' part of 7 ^ only. Collecting all presently 
available information, the anomalous dimensions for the resummed unpolarized singlet evolution 
in the DIS scheme are given by 
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Here 7o and 7l denote the leading and next-to-leading order singlet anomalous dimension ma- 
trices. 

The resummed terms beyond 0(a 2 ) in eq. (10) do not comply with the energy-momentum 
sum rule for the parton densities, which requires 

% q (N,a s ) + lgq (N,a s ) =0, <y qg (N } a s ) + -y 99 {N, a.) = . (11) 
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1.000 E+00 
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Table 1: The expansion coefficients 6| and 6f for the small-x resummed anomalous dimension 
7 9g . The latter quantities are given for four active flavors. For comparison to previously employed 
estimates [|Tj also the ratios of these coefficients are shown. 

This relation is satisfied by the fixed-order anomalous dimensions order by order in a s . The 
method to restore the sum rule with the least impact on the small-x results is a (diagonal) 
subtraction at iV = 2, i.e., the addition of appropriate 5(1 — x) contributions to the higher- 
order quark-quark and gluon-gluon splitting functions. We will label this prescription as (A) 
below. Other possibilities are the inclusion of somewhat less singular 1/(N — 1) terms, later to 
be superseded by explicit calculations. In this manner a rough estimate can be obtained of the 
possible effect of subleading small- x contributions to the higher-order anomalous dimensions, 
particularly in those cases where only the first term is known currently. We will illustrate this 
procedure by the prescription (D) given by 

7^(iV)^ 7jJ (iV)(l-2[iV-l] + [iV-l] 3 ) for ij = qq,qg,gq. (12) 

The terms oc [N— 1] in 7 9g are taken from eq. (|9|), hence only the terms oc [iV— l] 3 in this quantity 
are adjusted according to eq. (|Tl"l). 

We are now ready to discuss the numerical effects of the small-x resummations. For definite- 
ness, we choose the initial distributions of the MRS (A') global fit Both the gluon and the 



sea quark densities behave as x~ 0A7 for a;^0 at the starting scale of Ql = 4 GeV 2 . The evolu- 
tion is performed for four massless flavors, and also A^ 4 -* = 231 MeV is adopted from the MRS 
analysis. We stress that these results are mainly theoretical illustrations. Detailed data analyses 
would require some flexibility of the input gluon density at small x, which is only indirectly 
constrained by structure function data, as well as the inclusion of heavy-flavor mass effects. 

Figure 1 displays the evolution of the proton singlet quark and gluon (momentum) distribu- 
tions, xS(x, Q 2 ) and xg(x,Q 2 ), in the DIS scheme. Different resummation approximations are 
compared with the NLO results (full lines). The results with the new b%'^ terms omitted are 
marked by NLx„. 
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In the quark sector the resummation corrections using prescription (A) are very large, e.g., 
they exceed a factor of four at x — 1CT 5 and Q 2 = 100 GeV 2 . This huge correction is entirely 
dominated by the quarkonic (upper row) anomalous dimensions. Omitting the contribu- 
tions, and even ignoring the Lx terms (that case is not shown in the figure), has an impact of 
less than about 20% on the singlet distribution 0, ||. The quark evolution is, however, very 
much affected by possible subleading contributions to 7 g9 and j gq , even by such terms which 
are smaller than those now found for 7 9g . This is illustrated by prescription (D), where the less 
singular pieces actually overcompensate the effect of the leading I = 1 contribution. 

The effect of the new contribution 7W to the anomalous dimension is, on the other hand, very 
substantial for the gluon density: the results even fall noticeably below the NLO evolution, and 
also below our previous lower estimate (D) at NLx q accuracy ||. There is no convergence so far, 
but the analysis of the known LO and NLO results supports some hope that the inclusion of also 



two further subleading terms (series), / = 2 and / = 3, may lead to a sufficiently stable result |~3 



At this point the question arises how large a correction the presently unknown resummed 7^ 
term may introduce. Experience in NLO and the Lx and NLa^ resummations suggests that the 
impact of this contribution is of the order of 10% or less, for the evolution of both xS and xg. 

We now turn to the proton structure functions. As we are working in the DIS scheme, 
iS(i, Q 2 ) already reflects F 2 up to the non-singlet pieces which are not relevant at small x []. 
Hence we directly turn to Fl(x,Q 2 ). The resummed results, employing the parton densities 
shown in Figure 1 and the / = 1 resummation of the coefficient functions ||, 



C L (N) = a s Cl(N) + a 2 s Cl(N) ms + J2cL <k a s ^ 1 ^j , (13) 



are depicted in Figure 2(a). Here C® and Cl stand for the leading and next-to-leading order 
coefficient functions. Note that a full NLx calculation of Fl(x, Q 2 ) requires the knowledge of the 
presently unknown next-order resummed coefficient function even in the DIS scheme. 

The resummation corrections are exceedingly large at the lower Q 2 values shown. Here 
they are entirely dominated by the resummed coefficient functions. At very small values of x the 
resummation does even violate the condition Fl < F 2 at Q 2 ~ 4 GeV 2 , thus requiring more terms 
in the resummation or an adjustment of the input densities at Qq. At high Q 2 > lOOGeV 2 , due to 
the decrease of a s and the parton evolution, the effects of the quarkonic anomalous dimensions 
and the coefficient functions are of the same order. However, besides the anomalous dimensions 
also the coefficient function will receive subleading corrections, which are presently unknown. To 
estimate the possible consequences of these terms in Cx mentioned above, Figure 2(b) repeats 
the calculation illustrated by Figure 2(a), but with an estimate for those unknown contributions 
Cl — > Cx(l — 2[N — 1]) beyond next-to-leading order. As in the quark evolution, already a 
moderate correction can lead to an even drastic overcompensation of the I = 1 effect, calling for 
the evaluation of the next resummation contributions to Cl- 



Let us summarize: Recent results in refs. [BL Ol for the first time allow a determination of a 



subleading small- 2 resummed anomalous dimension, namely the 'irreducible' part of 7 59 in the 
unpolarized case. We have extracted the corresponding expansion coefficients in the usual DIS 
scheme. The coefficients of the subleading small-re poles turn out to be mainly of opposite sign 
than the Lx pieces, and they are typically much larger in iV-space. The numerical impact of 
these additional terms on the proton's parton densities and structure functions has been studied. 



lr The resummation of the small- a; terms (Lx) for the non-singlet structure functions was performed in | fj"7j . 
The corresponding corrections are smaller than 1% over the whole x range. 
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It is largest for the gluon evolution, where a substantial overcompensation of the positive leading 
resummation effect takes place, but less important for xE and Fl, which are dominated by the 
quarkonic anomalous dimensions and the coefficient functions. Note, however, that the 'energy- 
dependent' contributions still need to be derived. All in all, more terms need to be calculated in 
the small- x expansions, both of the anomalous dimensions and the coefficient functions, in order 
to arrive at stable resummation predictions. 
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Fig. 1: The small-a: evolution of the singlet quark and gluon densities including the resummed NLx q 
kernels [||] and the new terms of the gluon-gluon anomalous dimension |l2|] as compared to the 
NLO results. Two prescriptions for implementing the momentum sum rule have been applied. 
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F i g . 2: The small-x behavior of the longitudinal structure function for the parton distributions shown 
in Fig. 1. The resummed coefficient functions of ref. have been employed, in (a) without, in (b) 
with a moderate subleading contribution. 
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